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Abstract 

In the first stage of this work, we perform detailed calculations for the cross sections of the 
electron capture on nuclei under laboratory conditions. Towards this aim we exploit the advantages 
of a refined version of the proton-neutron quasi-particle random-phase approximation (pn-QRPA) 
and carry out state-by-state evaluations of the rates of exclusive processes that lead to any of the 
accessible transitions within the chosen model space. In the second stage of our present study, 
we translate the above mentioned e ~-capture cross sections to the stellar environment ones by 
inserting the temperature dependence through a Maxwell-Boltzmann distribution describing the 
stellar electron gas. As a concrete nuclear target we use the 66 Zn isotope, which belongs to the 
iron group nuclei and plays prominent role in stellar nucleosynthesis at core collapse supernovae 
environment. 
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I. INTRODUCTION 


Weak interaction processes occuring in the presence of nuclei under stellar conditions 
play crucial role in the late stages of the evolution of massive stars and in the presupernova 
stellar collapse jl-4]. As it is known, the core of a massive star, at the end of its hydrostatic 
burning is stabilized by electron degenerecy pressure as long as its mass does not exceed an 
appropriate mass (the Chandrasekhar mass limit, Mch) jd-S]. When the core mass exceeds 
Mch, electron degeneracy pressure cannot longer stabilize the center of the star and the 
collapse starts. In the early stage of collapse electrons are captured by nuclei in the iron 
group region 0 . 01 - 

During the pre-supernova evolution of core collapse supenova, the Fermi energy (or equiv¬ 
alently the chemical potential) of the degenerate electron gas is sufficiently large to overcome 
the threshold energy E t h r ( E t h r is given by negative Q values of the reactions involved in 
the interior of the stars) fj and the nuclear matter in the stellar core is neutronized. This 
high Fermi energy of the degenerate electron gas leads to enormous e _ -capture on nuclei 
and reduces the electron to baryon ratio Y e [10|, 111] • In this way, the electron pressure is 
reduced and the energy as well as the entropy drop. One of the important characteristics of 
the early pre-explosion evolution is the fact that electron capture on nuclei (specifically on 
nuclei of the pf shell) plays a key role jl2. 13|. 

In the early stage of collapse (for densities lower than a few 10 10 gr cm -3 ), the electron 
chemical potential is of the same order of magnitude as the nuclear Q value, and the e~ 
capture cross-sections are sensitive to the details of GT strength distributions in daughter 
nuclei. For this reason, some authors restrict the calculations only to the GT strength and 
evaluate e _ -capture rates on the basis of the GT transitions (at these densities, electrons 
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Various methods, 


are captured mostly on nuclei with mass number A < 60) 
used for calculating e~-capture on nuclei during the collapse phase, have shown that this 
process produces neutrinos with rather low energies in contrast to the inelastic neutrino- 
nucleus reactions occuring in supernova 


15 


18]. These neutrinos escape the star carrying 


away energy and entropy from the core which is an effective cooling mechanism of the 
exploding massive star Il9fl . For higher densities and tempratures, e~ capture occurs on 
heavier nuclei A > 65 |6, 8, Ill - |l3 ]. As a consequence, the nuclear composition is shifted 
to more neutron-rich and heavier nuclei (including those with N > 40) which dominate the 
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matter composition for densities larger than about 10 w gr cm 
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, 5,Q, 


20 ]. 


The first calculations of stellar electron capture rates for iron group nuclei have been 
performed by employing the Independent Particle Model (IPM) J2J]. Recently, similar studies 


have been addressed by using Continuum RPA (CRPA) 21|, large scale shell model 22, 23], 


RPA 9], etc 24], In the present work e _ -capture cross sections are obtained within a 
refined version of the Quasi-Particle Random Phase Approximation (QRPA) which is reliable 
for constructing all the accessible final (excited) states of the daughter nuclei in the iron 


group region of the periodic table 25-134]. For the description of the required correlated 


nuclear ground s; 
the BCS theory 


fates we c 
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n 


etermine single-particle occupation numbers calculated within 
|. Our nuclear method is tested through the reproducibility of 


experimental muon capture rates relying on detailed calculations of exclusive, partial and 
total muon capture rates 354401]. The agreement with experimental data provided us with 


high confidence level of our method and we continued with the calculations of electron 
capture cross sections in supernova conditions (where the densities and temperatures are 
high) using the pn-QRPA method. I 11 this paper, we performed calculations for m Zn isotope 
(it belongs to the iron group nuclei) that plays prominent role in core collapse supernovae 
stellar nucleosynthesis |l6{, [17|, 41], 

Our strategy in this work is, at first to perform extensive calculations of the transition 
rates for all the above mentioned nuclear processes assuming laboratory conditions, and then 
to translate these rates to the corresponding quantities within stellar environment through 


the use of an appropriate convolution procedure Qfl 0 Ir] . Gj |. To this purpose, we 
assume that leptons under such conditions follow Maxwell-Boltzmann energy destribution 

fl Gi- 


ii. CONSTRUCTION OF NUCLEAR GROUND AND EXCITED STATES 

Electrons of energy E e are captured by nuclei interacting weakly with them via W~ boson 
exchange as 

(A, Z) + e -A (A, Z — 1)* + v e (1) 

The outgoing v e neutrino carries energy E v while the daughter nucleus (A, Z — 1) absorbs 
a part of the incident electron energy given by the difference between the initial Ei and the 
final Ef nuclear energies as E u = Ef — E t . 
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The nuclear calculations for the cross sections of the reaction (P) start by writing down 
the weak interaction Hamiltonian 7 -L w which is given as a product of the leptonic, j l ^ pt , and 
the hadronic, currents (current-current interaction Hamiltonian) as 


n w = j" 


( 2 ) 


where G = Gfcos9 c with Gf and 9 C beingjthe well known weak interaction coupling constant 


and the Cabbibo angle, respectively 


[27, 28, 


42]. 


From the nuclear theory point of view, the main task is to calculate the cross sections of 
the reaction (P) which are based on the evaluation of the nuclear transition matrix elements 
between the initial |i) and a final |/) nuclear states of the form 


(f\H w \i) = -j=rJ d 3 xe~ i ^{f\JAi). 


( 3 ) 


The quantity stands for the leptonic matrix element written in coordinate space with 

q being the 3-momentum transfer. For the calculation of these transition matrix elements 
one may take advantage of the Donnelly-Walecka multipole decomposition which leads to a 
set of eight independent irreducible tensor multipole operators containing polar-vector and 
axial-vector components 42] (see Appendix [A]) . 

In the present work, in Eq. ([3]) the ground state of the parent nucleus |i) is computed 
by solving the relevant BCS equations which give us the quasi-particle energies and the 
amplitudes V and U that determine the probability for each single particle level to be 
occupied or unoccupied, respectively [ 27 ]]. Towards this aim, at first, we consider a Coulomb 


corrected Woods-Saxon potentia. 


of the strong nuclear field 45 


with a spin orbit part as a mean field for the description 


46]. For the latter potential we adopt the parametrization of 


IOWA group 47]. Then, we use as pairing interaction the monopole part of the Bonn C-D 
one meson exchange potential. The renormalization of this interaction to fit in the 66 Zn 
isotope, is achieved through the two pairing parameters g p '^ ir , p (n) for proton (neutron) 
pairs, the values of which are tabulated in Table HI 

As it is well known, the pairing parameters g^ r , are determined through the reproduction 
of the energy gaps, A^, from neighboring nuclei as (3-point formula) 

1 


A exp = 
n 4 


S n (A -1,Z)- 2S n (A, Z) + S n (A + 1, 2) 


( 4 ) 
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TABLE I. Parameters for the renormalization of the interaction of proton pairs, g p air , and neutron 
pairs, g™ air ■ They have been fixed in such a way that the corresponding experimental gaps, A p Xp 
and A(f p of 66 Zn isotope, are quite accurately reproduced. 


Nucleus 

Q n ■ 
Vpair 

V 

Qpair 

A n p {MeV) 

A the° r ( Me y) 

A p xp ( MeV ) 

Af eor {MeV) 

66 Zn 

1.0059 

0.9271 

1.7715 

1.7716 

1.2815 

1.2814 


A exp 
P 



1,Z-1)~ 2 Sp(A, Z ) + S P (A + 1 , Z + 1 ) 


(5) 


where S p and S n are the experimental separation energies for protons and neutrons, respec¬ 
tively, of the target nucleus (A,Z) and the neighboring nuclei (Ail, Z± 1) and (Ail, Z). For 
the readers convenience in Table nn we show the values of experimental separation energies 
for the target Zn and the neighboring nuclei %Cu, Ga, Zn and %Zn. 


TABLE II. The experimental separation energies in MeV for protons and neutrons of the target 
(A,Z) and neighboring (A ± 1, Z ± 1) and (A A 1, Z ) nuclei. 


Nucleus 

S n {A-l,Z) 

Sn{A,Z) 

S n {A + l,Z) 

S P {A — 1,Z — 1) 

S P (A,Z) 

Sp{A + 1,^ + 1) 

66 Zn 

7.979 

11.059 

7.052 

7.454 

8.924 

5.269 


Subsequently, the excited states | /}_ of the studied daughter nucleus m Cu are constructed 
by solving the pn-QRPA equations 25-34], which in matrix form are written as 25] 



= W 



( 6 ) 


42 jTr denotes the excitation energy of the QRPA state | JJ) with spin J and parity tt. 

The solution of Eqs. (ED is an eigenvalue problem which provides the amplitudes for 
forward and backward scattering X and Y, respectively, as well as the QRPA excitation 
energies 42 


27 


30]. In our method the solution of the QRPA equations is carried out 


separately for each multipole set of states | J n ). 


5 

































For the renormalization of the residual 2-body interaction (Bonn C-D potential), the 
strength parameters for the particle-particle (g pp ) and particle-hole (g p h) interaction entering 
the QRPA matrices A and B, are determined (separately for each multipolarity) from the 
reproducibility of the low-lying experimental energy spectrum of the final nucleus. The 
values of these parameters in the case of the spectrum of m Cu are listed in Table ELD 

TABLE III. Strength parameters for the particle-particle (g pp ) and particle-hole (g p h) interaction 
for various multipolarities (for the rest of multipolarities, J n < 5^, the bare 2-body interaction 
has been used), in the case of the spectrum of m Cu nucleus. 


Positive Parity States 

Negative Parity States 


0+ 1+ 2+ 3+ 4+ 

1" 2” 3' 4" 

9pp 

0.827 0.547 0.686 0.854 1.300 

0.994 0.200 0.486 0.622 

9ph 

0.336 0.200 1.079 0.235 0.200 

1.200 0.200 1.200 1.200 


At this point, it is worth mentioning that for measuring the excitation energies of the 
daughter nucleus m Cu from the ground state of the initial one 66 Zn, it is necessary a shifting 
of the entire set of QRPA eigenvalues. Such a shifting is required whenever in the pn-QRPA 
a BCS ground state is used, a treatment adopted by other groups previously 0, 0, m Q- 


The shifting for the spectrum of the daughter nucleus 66 Cu, is done in such way that the 
first calculated value of each multipole state of m Cu (i.e. 1+, 2\. ..etc.), to approach as close 
as possible the corresponding lowest experimental multipole excitation. Table IIVI shows 
the shifting applied to our QRPA spectrum for each multipolarity of the parent nucleus 
m Zn. We note that, a similar treatment is required in pn-QRPA calculations performed 
for double-beta decay studies where the excitations derived for the intermediate odd-odd 
nucleus (intermediate states) through p-n and n-p reactions from the neighboring nuclei, 


left or right nuclear isotope, do not match to each other 48, [49]. The resulting low-energy 
spectrum after using the parameters of Tables [T] and IIIII and the shifting shown in Table IIVI 
agrees well with the experimental one (see Fig. [T|) . We must also mention that, usually, 
in nuclear structure calculations we test a nuclear method in two phases: First through 
the construction of the excitation spectrum as discussed before, and second through the 
calculations of electron scattering cross sections or muon capture rates. Following the above 
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TABLE IV. The shift (in MeV) applied on the spectrum (seperately of each multipole set of states) 
of m Cu isotope, daughter nucleus of the electron capture on 66 Zn. 


Positive Parity States 

Negative Parity States 

0+ 

0.90 

0" 

5.00 

1+ 

2.50 

1“ 

6.80 

2+ 

2.55 

2" 

3.85 

3+ 

2.50 

3~ 

2.60 

4+ 

1.75 

4" 

3.55 

5+ 

0.55 

5“ 

3.00 


steps, we test the reproducibility of the relevant experimental data for many 
employed in nuclear applications (nuclear structure and nuclear reactions) 
astrophysics Q0- 


nuclear models 
and in nuclear 


III. RESULTS AND DISCUSSION 


In this work we perform detailed cross section calculations for the electron capture on 
m Zn isotope on the basis of the pn-QRPA method. The required nuclear matrix elements 
between the initial |./,} and the final \Jf) states are determined by solving the BCS equations 


for the ground state [251,1271, (28] and the pn-QRPA equations for the excited states 271430] 
(see Sect. ED- For the calculations of the original cross sections, a quenched value of qa 
(see Appendix [B]) is considered which subsequently modifies all relevant multipole matrix 


elements 


39 


M 


50 


5l| 


At this point of the present work and in order to increase 


method, we perform total muon capture rates calculations 35-4401]. The comparison with 


he confidence level of our 


experimental and other theoretical results is shown in Sect. IIII A1 Afterwards, we study 
in detail the electron capture process as follows: i) Initially we consider laboratory condi¬ 
tions, i.e. the initial (parent) nucleus is considered in the ground state and no temperature 
dependence is assumed (see Sect. IIII B ID . ii) Second, we consider stellar conditions, i.e. 
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FIG. 1. Comparison of the theoretical excitation spectrum (resulting from the solution of the 
QRPA eigenvalue problem) with the low-lying (up to about 3 MeV) experimental one for m Cu 
nucleus. As can be seen, the agreement is very good below 1 MeV but for higher excitation energies 
it becomes moderate. 

the parent nucleus is assumed to be in any initial excited state and due to the e _ -capture 
process it goes to any final excited state of the daughter nucleus. At these conditions it is 
necessary to take into account the temperature dependence of the cross sections (see Sec. 

mm [si. 

A. Calculations of Muon Capture Rates for 6b Zn 

Despite the fact that the muon capture on nuclei does not play a crucial role in stellar- 
nucleosynthesis, it is, however, important to start our study from this process since, the 
nuclear matrix elements required for an accurate description of the ^-capture are the same 
for all semi-leptonic charge-changing weak interaction processes. In addition, the excitation 
spectrum of the daughter ( A , Z — 1) nucleus, as we saw before, is in good agreement with 
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the experimental data. 

The calculations of the muon capture rates are performed in three steps: In the first 
step we carry out realistic state-by-state calculations of exclusive Ordinary Muon Capture 
(OMC) rates in 66 Zn isotope for all multipolarities with J 77 < 5 1 * 1 (higher multipolarities 
contribute negligibly). The appropriate expression for the exclusive muon capture rates is 
written as: 


A s «j F = Ajj = 2 G 2 {<S> u ) 2 R,q 2 f - £,,)||0+)| 2 + - -75"“"")||0+)| 2 (7) 


where $i s represents the muon wave function in the Is muonic orbit. The operators in 
Eq. ((71) refer to as Coulomb M.j, longitudinal Cj, transverse electric Tf and transverse 
magnetic / Jj iag multipole operators (see Appendix [A]) . The factor Rf in Eq. ((7j) takes into 
consideration the nuclear recoil which is written as Rf = ^1 + qj/M tar ^j , with M targ being 
the mass of the target (parent) nucleus. 

Due to the fact that there are no available data in the literature for exclusive muon 
capture rates, the test of our method is realized by comparing partial and total muon 
capture rates with experimental data and other theoretical results 39, |40]. Towards this 


purpose, our second step includes calculations of the partial /i~-capture rates for various 
low-spin multipolarities, A jn (for J 77 < 4A), in the studied nucleus. These partial rates are 
found by summing over the contibutions of all the individual multipole states of the studied 
multipolarity as 


A jvr = ^2 = 2 G' 


; ($ ls ) 2 [ ^2q 2 f R f \(Jj\\{Mj - Cj)\\0+ s )\ 2 

f 

+Y J q 2 f R M j im i - fr 9 n )iio+>i 2 
/ 


( 8 ) 


(/ runs over all states of the multipolarity | J 77 )). We also estimate the percentage (portion) 
of their contribution into the total p-capture rate for the most important multipolarities. In 
Table fV! we tabulate the individual portions of the low-spin multipole transitions (J 77 = 4A). 
As can be seen, the contribution of the 1“ multipole transitions is the most important 
multipolarity exhausting about 44% of the total muon-capture rate. Such an important 
contribution was found in 16 O and 48 C*a isotopes studied in Ref. [37]. 

In the last step of testing our method, we evaluate total muon-capture rates for the m Zn 
isotope. These rates are obtained by summing over all partial multipole transition rates (up 
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TABLE V. The percentage of each multipolarity into the total muon-capture rate evaluated with 
our pn-QRPA method. 


Positive Parity Transitions 

Negative Parity Transitions 

r 

Portions (%) 

J n 

Portions (%) 

0+ 

8.22 

0" 

7.94 

1+ 

21.29 

1“ 

44.21 

2+ 

2.85 

2~ 

13.32 

3+ 

1.58 

3" 

0.34 

4+ 

0.01 

4~ 

0.23 


to = 4 : ± ) as 


A « = E = E E A - 


(9) 


For the sake of comparison, the above mentioned /^-capture calculations have been carried 


out using the quenched value Qa = 1.135 


39 


40]. The results are listed in Table lYll where 


we a 


and 


so include the experimental total rates as well as the theoretical ones of Refs. 39 ] 


40]. Moreover, in Table IVT1 we show the individual contribution into the total muon 


capture rate of the polar-vector (A^ t ), the axial-vector (A^ t ), and the overlap (A^ 0 f) parts. 
As can be seen, our results obtained with the quenched qa coupling constant are in very 
good agreement with the experimental total muon-capture rates (the deviations from the 
corresponding experimental rates are smaller than 7%). This agreement provides us with 
high confidence level for our method. 


B. Electron Capture Cross Section 

After acquiring a high confidence level for our nuclear method, we proceed with the 
main goal of the present study which concerns the calculations of the electron-capture cross 
sections. As mentioned before, this includes original (see Section [ill B ID and stellar electron 
capture investigations fSection llll B 2D . 
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TABLE VI. Individual contribution of Polar-vector, Axial-vector and Overlap parts into the to¬ 
tal muon-capture rate. The total muon capture rates obtained by using the pn-QRPA with the 


quenched value of gA = 1.135 for the medium-weight nucleus bb Zn, are compared with the available 
experimental data and with the theoretical rates of Ref. 


39] and Ref. 


40|j 


Total Muon-capture rates A to< (xl0 6 s Q 



Present pn-QRPA Calculations 

Experiment 

Other theoretical Methods 

Nucleus 

\V \A aVA 

ly tot tot ly tot 

A tot 

a To! 

A tt or [39] 

a tr m 

m Zn 

1.651 4.487 -0.204 

5.934 

5.809 

4.976 

5.809 


1. Original Electron Capture Cross Section on bb Zn isotope 


The original cross sections for the electron capture process in the m Zn isotope are ob¬ 
tained by using the pn-QRPA method considering all the accessible transitions of the final 
nucleus m Cu. In the Donnelly-Walecka formalism the expression for the differential cross 


section in electron capture by nuclei reads 


8 ] 


da ec _ G 2 f cos 2 6 c F(Z, E e ) 
27r (2 Ji + 1) 


dd 




J> 1 


x {[1 - + bsm 2 ®} [\{Jf\\rr 9 U )\ 2 + \(Jf\\ff\\Ji)\ 2 } 

~ [ {£i + £f) (l - acos $) - d\2Re(J f \\Ty ia9 \\Ji)(Jf\\TJ l \\Ji)*} 
+ ^ W(£'e, E U ){(1 + acos$)\(Jf\\Mj\\Ji)\ 2 

J> 0 

+ (1 + acos& — 26.sm 2< f>)|(J/||£./|| Jj)| 2 
- [-(1 + acosi) + d}2Re{J,\\Lj\\ Ji){J,||Aij||J,) - }} 


( 10 ) 


where F(Z, E e ) is the well known Fermi function 16]. The factor W ( E e , E u ) = E 2 /( 1 + E U /M T ) 
accounts for the nuclear recoil jg , M F is the mass of the target nucleus and the parameters 
a, b, d are given e.g. in Ref. 27]. The nuclear transition matrix elements between the 


initial state | Jj) and a final state | J f) correspond to the Coulomb A4 jm , longitudinal Cjm, 
transverse electric Tf M and transverse magnetic Tj^ 9 multipole operators (discussed in 
Appendix [A]) 


11 
























From the energy conservation in the reaction (JT]), the energy of the outgoing neutrino E u 


is written as 


E v — E e — Q + Ei — Ef 


( 11 ) 


which includes the difference between the initial E % and the final Ef nuclear states. The Q 
value of the process is determined from the experimental masses of the parent (Mj) and the 
daughter ( Mf ) nuclei as Q = Mf — Mj Jt|. 

It is worth mentioning that for low momentum transfer, various authors use the ap¬ 
proximation q —> 0 for all multipole operators of Eq. ([TO . Then, the transitions of the 
Gamow-Tcller operator ( GT + = ^Vr+cq), provide the dominant contribution to the total 
cross section [7]. 

In performing detailed calculations for the original electron capture cross sections on m Zn 
isotope we assumed that (i) the initial state of the parent nucleus 66 Zn is the ground state 
|0 + ) and (ii) the nuclear system is under laboratory conditions (no temperature dependence 
of the cross sections is needed). The cross sections as a function of the incident electron 
energy E e are calculated with the use of realistic two-body interactions as mentioned be¬ 
fore. The obtained total original electron capture cross sections for m Zn target nucleus 
are illustrated in Fig. [2] where the individual contributions of various multipole channels 
(J n < 5 1 * 1 ) are also shown. The electron capture cross sections in Fig. El exhibit a sharp 
increase by several orders of magnitude within the first few MeV above energy-threshold, 
and this reflects the GT + strength distribution. For electron energy E e > lOMeV the calcu¬ 
lated cross sections show a moderate increase. From experimental and astrophysical point 
of view, the important range of the incident electron energy E e is up to 30 MeV. At these 
energies the 1 + multipolarity has the largest contribution to the total electron capture cross 
sections [7, s|. In the present work we have extended the range of E e up to 50 MeV since at 
higher energies (around 40 MeV) the contribution of other multipolarities like 1 _ , 0 + and 
0~ become noticeable and can not be omitted (see Fig. E]). 

From the study of the original electron capture cross sections we conclude that, the 
total cross sections can be well a ppr oximated with the Gamow-Tcller transitions only in the 
region of low energies [ 7 -10, 13, 14], For higher incident electron energies the inclusion of 
the contributions originating from other multipolarities leads to better agreement 
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FIG. 2. Original total cross sections of electron-capture on the 66 Zn (parent) nucleus calculated 
with pn-QRPA method as a function of the incident electron energy E e . The individual contribu¬ 
tions of various multipole channels (for J n < 5 1 * 1 ) are also demonstrated. 

2. Stellar Electron Capture on m Zn isotope 


As it is well known, electron capture process plays a crucial role in late stages of evolution 
of a massive star, in presupernova and in supernova phases [l 4 ]. In presupernova collapse, 
i.e. at densities p < 10 10 g cm~ 3 and temperatures 300 keV < T < 800 keV, electrons are 


10 


13 


captured by nuclei with A < 60 [7- 
p > 10 10 gcm~ 3 and temperatures T ~ 1.0 MeV, electron capture process is carried out on 


14]. During the collapse phase, at higher densities 


heavier and more neutron rich nuclei with Z < 40 and N > 40 


B 


11 


131- 


In an independent particle picture, the Garnow-Teller transitions (which is the most 
important in the electron capture cross section calculations) are forbidden for these nuclei 0 . 
However, as it has been demonstrated in several studies, GT transitions in these nuclei are 
unblocked by finite temperature excitations 19|, |20|. At high temperatures, T ~ 1.5 MeV, 


GT transitions are thermally unblocked as a result of the excitation of neutrons from the 
pf-shell into the g 9 / 2 orbital. 

For astrophysical environment, where the finite temperature and the matter density ef¬ 
fects can not be ignored (the initial nucleus is at finite temperature), in general, the initial 
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nuclear state needs to be a weighted sum over an appropriate energy distribution. Then, 
assuming Maxwell-Boltzmann distribution of the initial state |i) in Eq. (fIU|) [7, 0| , the total 
e _ -capture cross section is given by the expression 8] 

G 2 f cos 2 9 c (2 Ji + \)e~ Ei ^ kT ) 

a(E e ,T)= F o _ c J2 F ( Z ,E e y J 


2tt 

l 

x 'y ' (E e — Q + Ei — Ef ) 
f,J 


G(Z,A,T) 

moj\f )\ 2 

(2 Ji + 1 ) 


( 12 ) 


The sum over initial states in the latter equation denotes a thermal average of levels, with the 
corresponding partition function G(Z,A,T) js]. The finite temperature induces the thermal 
population of excited states in the parent nucleus. In the present work we assume that 
these excited states in the parent nucleus are all the possible states up to about 2.5 MeV. 
Calculations involving in addition other states lying at higher energies shows that they have 
no sizeable contribution to the total electron capture cross sections. As mentioned before, 


for the evaluation o 


gA is necessary 


39 


40|, 


re total electron capture cross sections, the use of a quenched value of 
51| . Since the form factor Fa{(^) multiplies the four components 


50 


of the axial-vector operator [see Eqs. (1A3|) - (1A6|) ]. a quenched value of gA must enter the 
multipole operators generating the pronounced excitations 0 _ , l ± ....etc. For this reason, in 
our QRPA calculations we multiplied the free nucleon coupling constant gA = 1.262 by the 
factor 0.8 HQQQ. 

The results coming out of the study of electron capture cross sections under stellar condi¬ 
tions are shown in Fig. [3]where the same picture as in the original cross section calculations, 
but now with larger contribution is observed. As discussed before, the dominant multipolar¬ 
ity is the 1 + , which contributes more than 40% into the total cross section. In the region of 
low energies (up to 30 MeV), the total e _ -capture cross section can be described by taking 
into account only the GT transitions, but at higher incident energies the contributions of 
other multipolarities become significant and can not be omitted. 

The percentage contributions of various multipolarities (with J 7r < 5 1 * 1 ) into the total 
e _ -capture cross section at T=0.5 MeV and for incident electron energy E e = 25 MeV are 
tabulated in Table IVUI In addition, at this Table we list the values of the individual e~ - 
capture cross sections of each multipolarity with < 5 ± . More specifically, for E e = 
25 MeV the 1 + multipolarity contributes about 44%, the 0 + contributes about 26% and 
the 1“ about 11%. The contributions coming from other multipolarities are less important 
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FIG. 3. Electron-capture cross sections for the 66 Zn parent nucleus at high temperature (T=0.5 
MeV) in stellar environment obtained assuming Maxwell-Boltzmann statistics for the incident 
electrons. The total cross section and the dominant individual multipole channels (J 71- < 5 1 * 1 ) are 
demonstrated as functions of the incident electron energy E e . 


(smaller than 5%). 

In performing state-by-state calculations for the electron capture cross sections, onr code 
has the possibility to provide separately the contribution of the polar-vector, the axial-vector 
and the overlap parts induced by the corresponding components of the electron capture 
operators. In Fig. [4] we illustrate the stellar differential cross sections of each individual 
transition of the polar-vector and axial-vector components. 


As mentioned before, our code gives separately the partial e _ -capture cross sections of 
each multipolarity. In order to study the dependence of the differential cross sections on the 
excitation energy lo throughout the entire pn-QRPA spectrum of the daughter nucleus, a 
rearrangement of all possible excitations oj in ascending order, with the corresponding cross 
sections, is required. This was performed by using a special code appropriate for matrices 
29]. In the model space chosen for m Zn isotope, for all multipolarities up to J 77 = 5^ we 
have a number of 447 final states. The differential electron capture cross sections illustrated 
in Fig. [4] present some characteristic clearly pronounced peaks at various excitation energies 
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TABLE VII. Total e _ -capture cross sections (in 10 ~ A2 MeV~ l cm 2 ) for E e = 25 MeV. The per¬ 
centage of each multipolarity into the total e _ -capture cross section evaluated with our pn-QRPA 
code also tabulated here. 


Positive Parity Transitions 

Negative Parity Transitions 

V 

^(xlO- 42 ^) 

Portions (%) 


-e(xl0- 42 ^) 

Portions (%) 

0+ 

31.164 

25.96 

0" 

5.288 

4.41 

1+ 

52.779 

43.98 

1“ 

13.409 

11.14 

2+ 

6.921 

5.77 

2~ 

3.262 

2.72 

3+ 

5.499 

4.58 

3~ 

0.905 

0.75 

4+ 

0.244 

0.20 

4~ 

0.299 

0.25 

5+ 

0.208 

0.17 

5~ 

0.042 

0.04 


oj. These peaks correspond mainly to 0 + , 1 + and 2 + transitions. More specifically, in the 
m Cu daughter nucleus the maximum peak correspond to the 0^ QRPA transition at u = 
2.538 MeV and other characteristic peaks correspond to l^lg" and 1^ 0 transitions, located 
at energies u = 3.194MeV, u = 3.686 MeV and u> = 6.555 MeV, respectively (see Fig. ED- 
There are also shown other less important peaks in Fig. [Q 

Before closing, it should be mentioned that, the e-capture cross sections presented in this 
work, may be useful in estimating neutrino-spectra arising from e-capture on nuclei during 
supernova phase. The knowledge of //-spectra at every point and time in the core is quite 
relevant for simulations of the final collapse and explosion phase of a massive star. As it 
is known [l{J, in the collapse phase, neutrinos are mainly produced by e-capture on nuclei 
and on free protons. The energy spectra of the emerging neutrinos from both reactions are 
important ingredients in stellar modelling and stellar simulations HQ- 

Furthermore, in core collapse simulations one defines the reaction rate of electron capture 
on nuclei given by 

R h = ^2 Y i A/ (13) 

i 

where the sum runs over all nuclear isotopes present in the astrophysical environment {Y l 
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FIG. 4. Individual contributions of the Polar-Vector, (Ay), and Axial-Vector, (A^), components 
as well as total electron-capture rate as functions of the excitation energy ui ( 66 Zn is the parent 
nucleus). 


denotes the abundance of a given nuclear isotope and A* is the calculated electron capture 
rate for this isotope). The rates of Eq. (fT3]l must be known for a wide range of the 
parameters: T (temperature) and p (nuclear density) of the studied star. Thus, for the 
calculation of the quantity Y ■ A of a specific nuclear isotope one needs to know in addition 
to the nuclear composition Y the electron capture rates A calculated as we have shown 
in our present work. The rates of electron capture on various nuclear isotopes and the 
corresponding emitted neutrino spectra in the range of the parameters (T, p, Y e ) describing 


o + 


1 


E. = 25.0 MeV 


Polar Vector 


66 Zn(e",v e ) 66 Cu 

T=0.5MeV 
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the star until reaching equilibrium during the core collapse, are comprehensively studied 


in Ref. 


HE, 


23] for a great number of nuclear isotopes by using the large scale shell 


model. We are currently performing similar calculations for a set of isotopes by employing 
the present pn-QRPA method 52], 

Furthermore, the average neutrino energy, (E u ), of the neutrinos emitted by e-capture on 
nuclei can be obtained by dividing the neutrino-energy loss rate (defined by an expression 
similar to Eq. (I13p by replacing the rate A, with the energy loss rate Ej) with the reaction 
rate for e-capture on nuclei Rh . Assuming e.g. power-law energy distribution for the neutrino 
spectrum produced by the e-capture in supernova phase, the average neutrino-energy (E u ) 
determines a specific supernova-neutrino scenario. In addition, the neutrino emissivity is 
obtained by multiplying the electron capture rate at nuclear statistical equilibrium with 


the neutrino-spectra 


19 


22 


23]. Finally we note that, the rates for the inverse neutrino 


absorption process are also determined from the electron capture rates obtained as discussed 
in this section 19]. 


IV. SUMMARY AND CONCLUSIONS 


The electron capture on nuclei plays crucial role during the presupernova and collapse 
phase (in the late stage e _ -capture on free protons is also significant). It becomes increasingly 
possible as the density in the star’s center is enhanced and it is accompanied by an increase 
of the chemical potential (Fermi energy) of the degenerate electron gas. This process reduces 
the clectron-to-baryon ratio Y e of the matter composition. 

In this work, by using our numerical approach based on a refinement of the pn-QRPA 
that describes reliably all the semi-leptonic weak interaction processes in nuclei, we studied 
in detail the electron capture process on 66 Zn isotope and calculated original as well as 
stellar e _ -capture cross sections. We tested our nuclear model (the pn-QRPA) through 
the reproducibility of orbital muon capture rates for this isotope. The agreement with 
experimental data and other reliable theoretical results of partial and total p-capture rates 
as well as of the percentage contributions of various low-lying excitations is quite good which 
provides us with high confidence level for the obtained cross sections. 

Our future plans are to extent the application of this method and make similar calcu¬ 


lations for other interesting nuclei 


52] - Also this method could be applied to other semi- 
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leptonic nuclear processes like beta-decay and charged-current neutrino-nucleus processes 
important in nuclear astrophysics and neutrino nucleosynthesis. 
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Appendix A: Nuclear Matrix Elements 

The eight different tensor multipole operators entering the above Equations (see Sec. 
IIIID refer to as Coulomb M.jm, longitudinal Cjm , transverse electric Tj l M and transverse 

j • r-rmaan 112 1 

magnetic / JM , are defined as 


MMqr) = Mffij + MS' 5 , C JM (qr) = l JM + P JU (Al) 

Tf M (qr)=ff M + Tf^, fZ , "(qr) = f']‘ff’' + fTM a ' a (A2) 

These multipole operators contain polar-vector as well as axial-vector parts and are writ¬ 
ten in terms of seven independent basic multipole operators as 


MS'(?r) = (q^)M(f(qr) 


(A3) 


Lifter) = ^MSV) 


(A4) 



g 

m n 




(AS) 
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jfmag( nr ) - 

L1 JM KV 1 ) ~ 


M, 


N 


Fr(rt)Ai I ( q r)-- f i v ( q l)Z' J M (qr) 


(A 6 ) 


iMjMiq r ) = jj 


N 


Fa(iDK(v) + i) + q 0 F P (qi))S"-’ M (q r 


2\\ 


(A7) 


iL JM(Q r ) = 


Fa{&) 


2 M l 


N 


Mi) 


S"m( 9 r) 


(AS) 


iTf M (qr) = F A (ql)J:' J M (qr) 


(A9) 


TTm 5 ^) = F A{ql)^ J M {qr) 


(A 10 ) 


where the form factors F x , X=1,A,P and /i v are functions of the 4-momentum transfer 
and Mtv is the nucleon mass. 

These multipole operators, due to the Conserved Vector Current (CVC) theory, are re¬ 


duced to seven new basic operators expressed in terms o: 


harmonics and vector spherical harmonics (see Refs. 27, (40|, |42j). The single particle re¬ 


spherical Bessel functions, spherical 


duced matrix elements of the form 0'i IIT/H. 7 ' 2 ), where Tf represents any of the seven basic 
multipole operators (Af^-, Q. J M , E^, S' A j, £"/, A^, A^) of Eq. flA3h - flA10jl . have been 
written in closed compact expressions as [27] 


{(ni(i)jil|r J |[(n 2 i 2 )j 2 ) = e-H ,« 2 V pfr 

li=0 


(All) 


where the coefficients P^ are given in Ref. [27|. I 11 the latter summation the upper index 
n max represents the maximun harmonic oscillator quanta included in the active model space 
chosen as n max = (Ah + N 2 — (3)/ 2, where N t = 2+ l u i=l,2, and f3 is related to the rank 


of the above operators 


27] 
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In the context of the pn-QRPA, the required reduced nuclear matrix elements between 
the initial |0+) and any final |/) state entering the rates of Eq. (J7J) are given by 

(/r-'lioj,) = £ 44EM [X hh ulvl + Y hh vl u”] (A12) 


J2>jl 


where Uj and v 3 are the probability amplitudes for the j-level to be unoccupied or occupied, 
respectively (see the text) 0.28, 

These matrix elements enter the description of various semi-leptonic weak interaction 


processes in the presence of nuclei 274341, 


3-44] 


Appendix B: Nuclear Form Factors 


In Eqs. (I A Ml) - flAlOJl the standard set of free nucleon form factors F x (ql), X= 1, A, P 


and /r (gjj reads 


Fi{<$ = 1-000 [l+( 


q 


1 -2 


'840 MeV 


(Bl) 


/(<) = 4.706 


l+( g ) 

v«/in 


F A (q;) = g A 


! + ( 


840 MeV' 


q ' 

1032 MeV' 


-2 


-2 


(B2) 

(B3) 


Fp(qt) = 


2 , _ 2Mn Fa^) 


(B4) 


q 1 2 3 + ml 

where Mn is the nucleon mass and gA is the axial vector free nucleon coupling constant (see 
the text). 
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